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Abstract 

Theories with curvature squared terms in the action are known to contain ghost modes in general. 
However, if we regard curvature squared terms as quantum corrections to the original theory, the 
emergence of ghosts may be simply due to the perturbative truncation of a full non-perturbative 
theory. If this is the case, there should be a way to live with ghosts. In this paper, we take the 
Euclidean path integral approach, in which ghost degrees of freedom can be, and are integrated 
out in the Euclideanized spacetime. We apply this procedure to Einstein gravity with a Weyl 
curvature squared correction in the inflationary background. We find that the amplitude of tensor 
perturbations is modified by a term of 0(a 2 H 2 ) where a 2 is a coupling constant in front of the 
Weyl squared term and H is the Hubble parameter during inflation. 



I. INTRODUCTION 

Previously in the literature corrections to the Einstein-Hilbert action formed with terms of 
higher order in the Riemann tensor have been considered. In the past these have been considered 
in connection with the possibility of a non-local theory |l( and also renormalisation . These 

higher derivative theories do, however, suffer from the problem of ghosts. There are, of course, 
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'good' ghosts, such as the Faddeev-Popov ghost; but higher derivative ghosts are 'bad' ghosts. 
Even if one adds 'small' higher derivative terms, so that the lagrangian is almost the same as the 
second order theory, the resulting theory is not well described perturbatively since there are now 
extra degrees of freedom. These bring with them a number of problems [l|, 

aaa. 

They are 

negative energy particles, with energies which blow-up in the limit where the higher derivative 
terms are removed from the Lagrangian. Ghosts also form states of negative norm and hence lead 
to non-unitarity. 

The simple harmonic oscillator can be made to mimic the effect of these terms by adding 
higher derivative terms to its Lagrangian. This was considered by Hawking and Hertog f|; they 
showed that despite these problems it is possible to obtain a sensible probability distribution for 
observations of the field. This result relied on the Euclidean path integral formulation of quantum 
theory. This, of course, echoes the No-boundary cosmology of Hartle and Hawking [8]. Other 
authors have considered the effect of extra terms in the lagrangian. In the cosmological case the 
full theory is not taken seriously, and instead back-substitution 

flQ 

is typically used. Here we 

run contrary to this trend by extending and comparing the result of Hawking and Hertog to the 
cosmological case where the Lagrangian contains a Weyl tensor squared term. 

It is normally the case for ghosts that even if one could deal with the non-interacting theory, the 
interacting theory would have the problem of run-away production of positive energy and negative 
energy particles. This would require further work, but is not expected to be a problem in the 
Euclidean approach where a late time boundary condition is specified jf]. [lo|. 

The very early Universe provides us with the ultimate high energy physics experiment. Post 
inflation the energies are sufficiently low that one could not rule out the occurence of higher order 
curvature terms in the Lagragian; however, during inflation these terms could be of importance. 
One could consider a whole series of terms of arbitrary order in the Riemann tensor, but it makes 
sense to start with just second order. 

Having accepted that there is the possibility of higher curvature terms, the question is how one 
should deal with them. One can handle them simply by back-substituting the modes from the 
usual equations of motion. This approach has been used in [5[|. Here one finds a momentum cut-off 
much higher than the Plank energy. 

On the other hand [11] argue that by using a perturbative approach, rather than solving the 
full system, one can miss essential physical features. For example, in the case of the heat equation 
initial data of compact support lead to a perturbative solution with the same support; however, 
the actual solution clearly spreads. 
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There is something special about the case where the only additional terms we consider are 
curvature squared terms. This theory has better renormalisation properties [jjj]. We deal with the 
slow roll case in four dimensions and therefore restrict our attention to the Weyl tensor squared 
terms. 

The structure of this paper is as follows. In $11] we review the simple harmonic oscillator. In 
§1111 we find the action. In §IVI we solve the eom. In $V]we find the wavefunction. In §VII we find 
the two point function; and in Will we draw our conclusions. 



II. THE HIGHER DERIVATIVE SIMPLE HARMONIC OSCILLATOR 

In [6| the authors consider a higher derivative simple harmonic oscillator, the extra terms being 
set up to emulate the situation in gravity. This has the advantage of fewer terms in the lagrangian 
and doing away with issues to do with gauge, thus providing clarity in the exposition of the method. 
One can expect that the method works in the case of gravity, and if one is not interested in the 
detailed form of the resulting observables it is more than adequet. However, as future experiments 
advance one may expect to be more interested in these details, which an ersatz theory cannot 
provide. Wishing to examine the form of the result more closely in gravity, the natural place to 
start is the tensor sector as it is already gauge invariant. 

Hawking and Hertog find that despite the problems presented by typical higher derivative 
theories it is possible to take the full theory seriously (not removing any modes) and arrive at a 
sensible probability distribution. Key to this method is the fact that they are considering the field 
as a Euclidean field. 

The action for a Euclidean higher derivative simple harmonic oscillator is 

I = fdr (^ 2 rr + i# + imV) ; (1) 

and, although this theory doesn't contain any interaction terms the authors argue that these would 
not cause further problems in the case where we are prescribing the field value on a late time surface. 
For real (p and real r this action is clearly positive semi-definite and hence results in a convergent 
path integral. The resulting equation of motion 

ot 2 <j>,TTTT — (f>,rr + rn 2 (f> = (2) 

has a general solution which is easily seen to be 

<f>(r) = A sinh(A+r) + B cosh(A+r) + C sinh(A_r) + D cosh(A_r) (3) 
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where 



A± = y 1 T \/l - 4m 2 a 2 

V2a 2 



m 

(4) 

a 



for small alpha (which is our area of interest since this results in a theory which has 'small' higher 
derivative terms). The ground state wave function at r = 0, ^ T =o, is the amplitude to go from zero 
4>(t = -co) = 0, </>, t (t = — oo) = to prescribed values 4>o, 0o,r at r = 0. The authors calculate 
this using the stationary phase approximation (and are thus assesing on a classical solution) , 

(/ 777 \ 2t77^ 777 ! OL \ 

-F [4>l T + -4) + — _ ^ <p <Po,r J (5) 

where 

F = 1 - 4m 2 « 2 

2a 2 (A+ + A_)(A_ - A+) 2 ' ( ' 

The Euclidean conjugate ground state wave function ^* is defined in a similar way, but integration 
goes from r = to oo; it is the analogue of complex conjugating the Lorentzian wavefunction. The 
associated probability is thus 

P(cP ,cP , T ) = = iV 2 exp (-2F (f^ + ^0 2 )) , (7) 

and we see that we may integrate over the unobserved 0o,r- If we had been working with cfioj, where 
t is the Lorentzian time, the sign infront of the (j)^ t term would have been wrong and integration 
not possible: the possibility is afforded by working in the Euclidean formalism. This results in a 
probability 



2Fm ( 2mF 2 

P(M = \ exp cp 2 (8) 

na \ a 



fm / 



ma\ . ,,. 

1 + — J exp(-m(l + ma)4>l) , (9) 



the approximation holding for small a. So we see that the second order theory is corrected by 
terms in a, rather than a 2 as we would see had we back-substituted the modes from the second 
order theory into the fourth order theory. 



III. PRELIMINARIES 



We are interested in actions with terms up to second order in the curvature tensor. Since we 
are working in four dimensions the Gauss-Bonnet term is a total derivative and so we may replace 
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R^R^ with a combination of R 2 and C a ^ u C a/3lMl/ . The R 2 term leads to an additional massive 
scalar field which doesn't cause any problems. 1 We are interested in the negative energy ghost 
component, which Boulware [lj| has shown to be contained in the tensor part of the Weyl squared 
term. Thus, we focus on this. 

The York-Gibbons-Hawking boundary term Bygh is usually added to the Einstein-Hilbert 



action in order to make t 



surface) well posed 



re variational problem (where one chooses the metric on a spacelike 



This comes down to the fact that if we are specifying 'q' at initial 
and final times we do not also wish to specify £ q.' In the context of the higher derivative theory 
we condsider here we will specify l q J and 'g,' so at first glance it might seem that we do not need 
Bygh', however, since this term also makes the Euclidean action for the inhomogeneous modes 
positive definite it transpires that we will need it. 

At zeroth order in slow roll the spacetime during inflation reduces to de Sitter space. We use 
the metric 

ds 2 = e 2p (-dr] 2 + (Sij + jij^dxi) , (10) 

where 7^ is a transverse traceless perturbation and we can conveniently take e 2p to be 1 or l/(Hr/) 2 
according to whether we wish to discuss flat or de Sitter space. From Boulware [jij we have 

^ ^aB^v = S(y MJfcOZ + 4 ^ <^0klm- (Hj 

This leads to 

C a ^C aP , v = \e- Ap {Wl 3 + 7wm)(-$ + lii,mm) ~ 2(V y , fc - 7^0(7^ - -/ifcj)) (12) 

1 

~ 2 



~< 4p Wijlij + 2 Hj,nnlij + Hjlij,kk + 7y,nn7ii,mm) , (13) 



at second order, where spatial boundary terms have been dropped. 2 Here and in what follows the 
repeated spatial indices i,j, k, . . . are to be summed over. Since the metric is conformally flat the 
Weyl tensor is vanishing at zeroth order and hence it is only necessary to calculate it to first order 
if one is finding the C 2 term to second order. 

Considerations so far lead us to start with the action 

S = f ^gd A x(^R-A- a 2 C a ^ v C a ^ u ) + B YGH . (14) 



1 One might expect a problem with the R 2 term; however, the only metric component having a non-degenerate 
higher derivative term due to this is the gravitational potential energy in the usual theory. The extra modes due 
to the higher derivatives therefore have positive energy. 

2 We are only ever concerned with the fields on future and past boundaries of a spacetime volume in the definition 
and propagation of the wavefunction. 
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where 



Bygh = J d 3 xVhK : 



(15) 



however, we will need an additional boundary term relating to the C 2 term. 

A. Action for fiat space 

If A = 0, the Minkowski spacetime is a solution of the theory. On this background the action 
([Til) becomes 



S 



dr/d 3 x 



^{lijlij ~ lij,klij,k) 



(16) 



at second order in the gravitational wave perturbation. We have dropped spatial boundary terms 
as these vanish for localised fields. We only consider the transverse traceless perturbations, and 
at second order these decouple from the scalar and vector perturbations. This gives us the ghost 
content of the theory. Considering now iS (the argument of the path integral defining the wave- 
function) rotated into the Euclidean with r = irj being the Euclidean 'time' we see that the only 
term preventing the Euclidean action / = — iS being positive definite on real Euclidean fields is 



(17) 

This leads us to the addition of a boundary term on the initial and final boundary surfaces. In 
this way we can do an integration by parts in the time coordinate and have a positive definite 
Euclidean action, so the path integral will be well defined. This boundary term is the same in the 
de Sitter case. 

B. Action for de Sitter space 



The action in the case of the metric (1101) is 



dr]d 3 x 



_ 3e 2p p '2 _ e 4p A 



a>- 



7ij7ij + 2 lijlij,kk + 7ij,nn7ij,mm 



(18) 
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where the background equation of motion (eom) holds. The eom for the tensor fluctuations is thus 
- {(e 2 %)' - e 2 P lijM ) + a 2 (7™ - 2rfl jM + 7ijl _) = . (19) 



4 

So expanding 7 as 

d 3 k 



^ = J 7j^£4 (^W^ (20) 



s=± 

where e s u (k) = = tfeffi), e*-(-fc) = ef/(£) and e£(fc)ef/(£) = 2S st we find 



Classically 7^ is taken as a real field. Then the second order action may be written in terms of 
the fourier components as 



s=± J 



'k 1 1 'fc 1 



(22) 



which we regard as an action for the real and imaginary parts of 7| as independent fields. 



C. Canonical formalism 

In (| it is argued that whilst most authors considering the canonical formalism for fourth order 
theories would take 'g' and l q > to be the canonical 'position' coordinates, in the path integral 
formalism one should describe a state by i q' and l q > at initial and final times. With this in mind 
we consider Q 7 | = 7^ and Q^ s i = 7^; these have conjugate momenta 

p i: < 24 > 

and similarly for the imaginary components of the fields. Thus we have a Hamiltonian 

K\ = 'All "| 2 + 2^( 7 f *7|') + ^|7|'| 2 + 2a 2 A; 2 |7f , | 2 + ^|7l| 2 + « 2 A ; 4 |7f| 2 (25) 

= K(P£Q - -^\P Y >\ 2 - Vl^'l 2 + ^\Q^\ 2 ~ 2a 2 fe 2 |Q f + a 2 fc 4 |Q 7S | 2 ,(26) 

'k in Aa z 7 fc 4 ~g 4 fc ~s fc 

where again we take the independent variables to be the real and imaginary components of each 
field. 

Taking the classical Hamiltonian (j26[) one sees the ghost instability is present because the only 
occurence of P 7 is in the PyQy term and hence the Hamiltonian can be made arbitrarily negative 
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by fixing Qy 7^ and taking P 7 appropriately large positive or negative. Unlike the case of a 
particle orbiting in a central potential, 3 this problem is present for a large volume of phase space, 



and it is this which makes the difference 



17(. Quantizing the higher derivative theory will not 



result in a lower bound on the energy of states. 



D. Flat space wavefunction 

In flat space the eom f|21|) becomes 

= U" + \ k ^i + ° 2 U" + 2k *< + *i) • < 27 > 

which is easy to solve in terms of exponentials. It factorizes to give solutions with k+ = k and 
k- = v^ 2 + l/(4a 2 ). 

The Lorentzian flat space wavefunction for the mode 7S is 

*|(Q 7 |,Q 7e =iV(77)exp | - a 2 k+k_(k + + fe_)|Q 7? | 2 



A- 



- za 2 A;_MQ; s Q + Q 7 §Q!sO + « 2 (^+ + MIQ 7 *'| 2 , (28) 

which may be obtained through the Euclidean path integral prescription, followed by rotation back 
to Minkowski time, and satisfies the Wheeler-de-Witt equation, 

5tf*| = (29) 

with 

<K\ l = MP*Q ,,)- JL|p s ,| 2 -I|Q J 2 

A; 2 

+— |Q 75 | 2 -2a 2 ^ 2 |Q ,,\ 2 + a 2 k 4 \Q- yi \ 2 , (30) 

where, as usual, P = —i8q. With the boundary term chosen to make the path integral defining 
the wavefunction well defined we will see that this is also true in the de-Sitter case. 



3 Of course, it is also the case that the classical hydrogen atom Hamiltonian can be made arbitrarily negative by 
taking the electron to have small momentum and be close enough to the nucleus. This problem with a small region 
of phase space suggests that the hydrogen atom has states of arbitrarily negative energy. Heisenberg's uncertainty 
principle stops a state being localised on this problematic region, and so there is a ground state. 
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IV. SOLVING THE EQUATION IN THE DE-SITTER CASE. 

The equation of motion (f2T1) for 7S with e 2p = H ] > 2 nas the two convenient factorizations, 

/ d 2 2 d 1 \s d 2 2 d \ , 

=(^2+-- + 1 -77^)(^2--- + 1 hf, (31) 



\dz 2 z dz Af3z 2 '\dz 2 z dz 

rid? 2d 1 \ / 9 d 2 d 9 1\„ . . 

" (f3? " ?S + ?) ( 2 a 5? " 22 S + 2 + 2 " ■ < 32 > 

where 2; = — krj and /3 = —H 2 a 2 . Hence the solutions vanishing in the upper half 77 plane are 

(1 + iz)e~ lz , (33) 

z3/2 ( J ivW z) -^ivW z) ) : ^ 

these each also solve the first factor in each of f)31j) . f)32j) . Whilst the mode (|33[) is conveniently 
the same as in the Einstein-Hilbert case, we will see that it differs in normalisation. One can also 
see that this mode is an eigenfunction of the first factor of (|32h with non-zero eigenvalue and thus 
obviously solves ([32]) . One can see from the second order equations these satisfy that they both 
also obay an equation of the form, 

?('!'" 4'*) < 35 » 

This, and similar equations, are useful in the normalization of the modes. 

A. Normalizing the modes 



With the canonical 'position' coordinates Q 7ij and Qy , the conjugate momenta are 

P % = -a 2 ^ . (37) 

The position space Hamiltonian these give rise to correctly generates the canonical Hamiltonian 
evolution equations, and thus these put the Poisson brackets in canonical form. The operator 
versions of these satisfy the equal time commutation relations 

[Q^ (77, x), P li4 fa, y)] = 2z5^ (x - y) , (38) 
[Qj Xv,x), P> . (v, f)] = 2«<5 (3) (x-y), (39) 
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where the factor of 2 occurs because summing over i, j sums over the "plus" and "cross" modes. 
If we use 

7 | = u k ag + u%a s _ % + v k bt* + v* k b s _ % (40) 

in (|20f) the normalized modes, 

u k = . H ^{1 - ikrj)e %kr) , (41) 



w fc = ; ^ ,/- P ^V /I +^(-^)lff (2) _(-M (42) 

Vfc 3 (l-8/3)V2 1 /j i^ITp ^ 

result in 

[al,a t J] = (2n) 3 5 st 5®(k-fy, (43) 

[6|,6| t ]=-(27r) 3 <5^ 3 )(fe-fc). (44) 

Thus the states created by 61 are of negative norm, as expected for ghosts. 



V. THE WAVEFUNCTION FOR A DE SITTER BACKGROUND 

In the case of de Sitter space the definition of the wavefunction in the Euclidean path integral 
formalism would normally be done in global coordinates with a foliation by equal time slices which 
are copies of S 3 . This is the Hartle-Hawking wavefunction [8]; however, there is a simplification in 
our case. We know [3] that for wavelengths passing through the horizon sufficiently late, 4 it suffices 
to use the coordinates of (|10p which cover only half the spacetime, resulting in a wavefunction 
similar to, and extending, that in [191 ] . So, the observables closely approximate those found in the 
calculation with S 3 hypersurfaces, and for all but the t < 20 on the S 2 of last scattering there is no 
loss in not doing the calculation in global coordinates. At late times the curvature terms don't affect 
the background equations of motion and we may use the coordinates with hypersurfaces which are 
copies of M 3 ; this is where we use the Euclidean formalism to integrate out over the unobserved 
variable Qy before rotating to Lorentzian time. The Wheeler de Witt equation, being based on 
a Hamiltonian formulation, is Lorentzian. It is generally convenient to give our expressions in a 
Lorentzian form and point out the changes in the Euclidean form where relevant. 

This wavefunction describing fluctuations about a de Sitter background may be found by ass- 
esing the action on a solution of the e.o.m. which has prescribed values of Q 7 , Qy (on a late time 



4 That is, sufficiently late so that they do not see the curvature of the universe at horizon exit. 
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surface at rjo) and vanishes in the upper half rj plane. This last condition restricts us to the modes 



in ||33D, ([MD, thus the 72 in PU) is 



7r 



Q-ys'(u(rj )v(r]) - v{rf )u(r])) - Q^Ju'^v (77) - v' (r] )u(r])) 



(45) 



fc u(rj )v' (rj ) - v(rj )u'(rj ) 

From this it is clear that any normalization of the modes will cancel out and thus doesn't concern us 
here. Reality of the field and its time derivative at 770 require Qy = (Q 7 s)* and Q 3 i = (Q s >)*- 
We are assessing the action on a solution of the eom and hence we are dealing with a boundary 
term on the future boundary: 5 



iS\ 



s=± 

where the coefficients 

A = 13 




S _i_ O ( ' 8 _.S 3 _.S j_ 0-,5_.S 

4z 2 ~k \Jk,zzz^-k ^k,zz^-k,z 'k'~k y 



\U,zzz^,z ^ \zzz^ ,z) 



UV z — VU 



so 



B 



fit + 2u?j 2 — 2?JU , + M,W 



V,zU,zz 



UV z — VU z 



20 



Wit 22 — It V • 



are obtained by the insertion of (1451) into (|46p . This leaves us with a wavefunction 



(46) 
(47) 



(48) 
(49) 
(50) 



k 3 
H 2 



N(t] ) exp 



—1 \q-?^-^(q^+q:^i\b+^\q_^c 



% 2 



The Wheeler-de-Witt equation, constructed from 
equivalent to 



1 B 2 - 

4^ + 4^- /? = ^ 



■Km = id^qt and P 



(51) 



-i8q, is then 

(52) 
(53) 
(54) 



which are satisfied by (@HD, (09]), (|5t)D . 



It is noted in [20( that the pre-exponential term in the wavefunction in the quadratic case is independent of the 
canonical coodinates; this is borne out by our result. 
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A. The small (3 approximation 

We will need to calculate A, B, C of (|4"8"]h (fl9l) . (|50l) . though little useful progress can be made 
employing the full form of the mode (I42j) . This is facilitated by only considering the physical case 
where (5 ~ (since H ~ and one can expect a is of order one). 

Once derivatives are neglected we are only dealing with expressions which are homogeneous of 
degree zero in u and v, so we can drop overall factors and use u, v in the form (j33j) . (j34|) . This can 
be further simplified by using 

R?) J. A (.-)-."»J >( z) 
A -zsm(A7r) 

where A = -y/l + 1//3/2, and dropping the factor of z/ sin(A7r). With real z and imaginary A the 
complex conjugate of J~\(z) is J\(z). Since /? is negative the e tnX suppresses the second term in 
(I55h more than any power of (3. The late time (small z) expansion for the Bessel functions 

Z X / z^ \ 

Jx{z) = 2*T(1 + A) V 1 " 4(1 + A) + 32(1 + A)(2 + A) + ' ' ' ) ' (56) 

shows us we can drop another factor, and we are left with the important part of the ghost mode 
being 

-XI - 



( 1 "4(r3Ay + 32(l-A )( 2-A ) + ---)' (5?) 



The first few terms of the resulting late time, small (5 series are 



i = -^ + i Q-2/3j + i(l + 2/3(-l + y / l + i) ) : + .... ,oSl 



2 



'-3 + J1 + 4 



'1+1-1)0 /3(3 + A /l + J 

VI. OBSERVABLES FROM TENSOR PERTURBATIONS 
A. The usual gravity 

Here we make a few notes about gravity where the Lagrangian consists only of the Einstein- 
Hilbert term. From equation (|18|) the second order part of the action is 



S 



dr]d 3 x^e 2p (jijTij - lij,klij,k) , (61) 
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with Fourier components as in equation (|22p given by 

%=E/^d<i 2 -^ 2 i7r)- (62) 

In this case, of course, we have canonical coordinates = and P~" = he 2p j s . ; and the 

'fcr fcr 'kr z kr 

Hamiltonian is given by 

*S = -rhl \ 2 + ~^hl\ 2 



k 4 1 'k 



e 2 Pk 2 



The wavefunction is constructed as in section (jVJ), where equation (|45p takes the simplified form, 

7§=^HA (64) 
and is plugged into a simplified form of equation (pt6j) , 

^ = -E / (^4l¥ 7 t 7 -^ ' (65) 

to give a wavefunction 



From this we see that the probability distribution for Q^a 

'k 

P(Q n ) = l^l 2 = l^exp { 2H ^f + z2) \Q^\ 2 ) (67) 

freezes out at late times (z 1 0); hence, observables involving Q 7 £ all freeze out. We note that 
in order for the various correlation functions of 7 to freeze out only the real component of the 
argument of the wavefunction (j66|) needs to freeze out. This is also observed in [19|]. We will see 
that a similar situation arises for the two point function when we include the Weyl squared term. 6 



B. With the Weyl squared term 

The probability distribution the wavefunction of equation (|5ip leads to is 

P(Q^Q^>) = |iV(r/o)| 2 exp(^ A ;|, + (z5 fc H r)0 )*) 



6 Linear combinations of the classical mode solutions in equations (|33|) . (|34|) freeze out since they start with a 
constant and have no term of order one in 2, thus we must have freeze-out of observables involving 7. 
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W(m)r exp ( — 2 ( |Q 7| |M - — ( q^,q; ? + q; £ ,q 7 , ) b + ^ |q 



|iV(r ?0 )| 2 exp 




4C 



Q, 



h 

kB 



1 



k 2 'I 



2 C 



2C 7fc 



(68) 



where ,4 = 23?(zA), etc. 

As before (section IV Aj) we may do a late time expansion of the ghost mode to select the 
dominant terms when (3 is small and negative, 



|JVfo)| 2 exp(|^( 



1 



1 



p z 



- + A(5+y-- I.;] :-' + 
10? 



\Q,t\ 2 



\/ 1 + ^ 2z | 2(3(1 + 10/3)z 2 



-1 + 3/3) 



(-1 + 



+ . 



+ . . . Qy> 



! + * 

Since — 1 < /3 < and z > it looks like the dependence on Qy is of the wrong sign to be 
integrated out. This problem does not emerge in the Euclidean formalism (as noted in [6j). We 
can rotate into the Euclidean and thus get an extra minus sign, which allows us to integrate over 
the unobserved variable. This leaves us with 



P(Q 7 ) = |iY(r ?0 )| 2 exp 



(69) 



2 V 2 

The simple fc dependence here is due to the fact that the k dependence can be removed from the 
eom that the modes satisfy; in the case of the flat space wavefunction the k dependence is not an 
overall factor. So observables freeze out at late times, with the two point function having a leading 
contribution for large k of 



H 



<hf>=^(l~SH^ + ...), 



(70) 



for the tensor perturbations 



VII. CONCLUSIONS 



We have shown that in higher derivative theories of gravity it is possible to live with ghosts, 
so that there is no need for back-substitution. So the result of taking the theory with the Weyl 
squared term seriously is not that we see corrections at order O(H), as we might have expected 
given the results for a higher derivative simple harmonic oscillator; instead, as in the case of back- 
substitution we get corrections at order 0(H 2 ). Thus, if we only consider the addition of a Weyl 
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squared term to the action we see no qualitative difference between back-substitution and taking 
the full theory seriously, as we do here. This similarity may not extend further: further work is 
required to determine whether the unusually high cut-off found by back-substitution [5J would be 
removed in our approach. 
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